arXiv:1508.05677v2 [math.PR] 22 Nov 2015 


NON-FIXATION FOR CONSERVATIVE STOCHASTIC DYNAMICS ON THE LINE 


RIDDHIPRATIM BASU, SHIRSHENDU GANGULY, AND CHRISTOPHER HOFFMAN 


Abstract. We consider Activated Random Walk (ARW), a model which generalizes the Stochastic Sandpile, 
one of the canonical examples of self organized criticality. Informally ARW is a particle system on Z with 
mass conservation. One starts with a mass density /r > 0 of initially active particles, each of which performs a 
symmetric random walk at rate one and falls asleep at rate A > 0. Sleepy particles become active on coming in 
contact with other active particles. We investigate the question of fixation/non-fixation of the process and show 
for small enough A the critical mass density for fixation is strictly less than one. Moreover, the critical density 
goes to zero as A tends to zero. This positively answers two open questions from [12,26]. 


1. The model description and main result 

Self organized criticality (SOC) is a universal property describing systems that have critical points as an 
attractor. In most of the examples of SOC simple local moves give rise to complex global properties. These 
systems are driven under their natural evolution to a state at the boundary between stable and unstable states 
without any fine-tuning of parameters. A canonical example of such a model with mass conservation is 
the Abelian Sandpile model which was proposed by Bak, Tang and Weisenfeld [3]. Additionally Dhar [7] 
proposed a collection of deterministic cellular automata as models exhibiting SOC. 

Several variations of the Abelian Sandpile model have been studied in the past with a view to better under¬ 
stand the critical behaviour. These variants have shown that adding randomness can lead to a qualitatively 
different critical behaviour. The best known example in this random setting is the Stochastic Sandpile 
model (SSM). On the graph Z this model starts with particles at every site distributed independently from a 
fixed distribution with mean /x. Any unstable site (a site with at least two particles) topples by emitting two 
particles which take independently one step of the simple random walk. The most important phenomenon 
in the study of SSM is the phase transition related to the long term fixation/non-fixation of the system. Em¬ 
pirically when the density is sufficiently small one observes activity in any finite window about the origin 
only finitely often, while when the density is sufficiently large the system does not fixate. 

Getting nontrivial bounds on the critical behavior has turned out to be extremely challenging. To make 
things mathematically tractable a closely related family of continuous time interacting particle systems 
known as Activated Random Walk was introduced. This system consists of particles which are in one 
of two states, active or sleepy. Activated random walk starts with particles distributed independently on the 
sites of a graph from a distribution with mean /i. Initially all of the particles are in the active state. As the 
system evolves each active particle performs a continuous time nearest neighbour symmetric random walk 
at rate 1 until it falls asleep, which happens with rate A. A particle in the sleepy state does not move but it 
becomes active whenever an active particle occupies its site. We denote this system on Z by ARW(/r, A). 

Note that the case when A is infinity is a variant of SSM. For a fixed sleep rate A, as the particle density /r 
increases, it is expected that the system shows a transition from almost sure local fixation to staying active 
forever almost surely. One of the first mathematically rigorous results about ARW was established in [26] 
where it is shown that for every A > 0 there is a critical particle density /Tc := /Uc(A) G 
ARW(^, A) locally fixates almost surely when // < ^c(A) and stays active almost surely when ^ > /ic(A). 
Quoting from [26, Section 7], 

“A proof that pc < 1/or the SSM and ARW remains as an open problem in any dimension. 

For the ARW, it should hold for all A, and moreover pcW —)• 0 A —)• 0. Yet, even a proof 

that /rc(A) < 1/or some X > 0 is missing.” 
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These problems have been re-iterated in [12], [6] and [32]. In this article we only look at ARW and for easy 
reference purpose we now record below the aforementioned conjectures in our setting. 

Conjecture 1. For ARW on Z, for all sleep rate A > 0, the critical mass density /ic(A) < 1. 

Conjecture 2. In the same set up as Conjecture 1, /Uc(A) —)• 0 as A —)• 0. 

The main result of this paper immediately provides a positive resolution of Conjecture 1 for small A. 

Theorem 1. Given /r > 0 there exists A^ > 0 such that ARW(/r, A) on Z with particle density fi and sleep 
rate A stays active almost surely for all A < A^. 

It was also established in [26] that //c(A) in a non-decreasing function of A. As a consequence we obtain 
the following corollary resolving Conjecture 2: 

Corollary 2. For ARW(/i, A) on Z, the critical particle density fid^) —)• 0 as A —)• 0. 

Remark 3. For any A > 0, it is easy to heuristically explain why ficW should be at most one. If // > 1, 
then on average “there are more particles than sites”, and since at most one particle can fall asleep at a site, 
the system should not fixate. In different settings this argument is formalized in [ 1 , 26, 30] , see Section 1 . 1 
for details. However, establishing /ic(A) < 1 requires understanding how activity is being sustained forever 
due to particle interaction, i.e. sleepy particles being woken up by active particles over and over. This makes 
the analysis substantially more difficult. 

1.1. Background. In the recent years, studies in non-equilibrium statistical mechanics have offered up 
a number of mathematically challenging interacting particle processes which exhibit phase transitions far 
away from equilibrium. A particular class of conservative models that has drawn significant attention is one 
where even though the mass is conserved, particles can exist in one of the two states: active and inactive. 
Typically an inactive particle becoming active requires interaction with one or more active ones. Paradigm 
examples in this class includes conserved lattice gases [29] and sandpile models with stochastic update rules 
[10,22,23]. In infinite volume, these models are believed to exhibit absorbing-state phase transition from an 
active phase as some model parameter (typically particle density) is varied. In finite volume, when run with 
a carefully controlled driven-dissipative mechanism, these systems are believed to exhibit the phenomenon 
of self-organized criticality [9, 11], where the system is attracted to a critical state, even though it is not 
explicitly tuned to this critical value. 

The transitions in these models are believed to belong to a universality class, referred to as the Manna 
class. Whether the Manna class exists as an autonomous universality class separate from the universality 
class of directed percolation (DP) seems not to have a broadly accepted answer among physicists at this 
point [4,19], it appears beyond the state-of-the-art techniques to obtain mathematically rigorous results on 
the critical or near-critical behaviour of these models (see [6] for some progress). Even the more basic 
questions about existence of phase transitions seem challenging and has only been settled in a few particular 
cases. One of the main challenges in studying these systems is the intricate long-range interaction caused 
by the conservation of particles, which makes it harder to apply some of the standard techniques in rigorous 
statistical mechanics. 

Activated Random Walk (ARW) is one of the well-known example in the above general class of models. 
As mentioned at the beginning, a major motivation to study ARW comes from the following Stochastic 
Sandpile Model (SSM), a variant of Manna’s model [22,23]. In the SSM on the line, started with an initial 
particle configuration of product measure with density p, a site with an isolated particle instantaneously 
becomes inactive, whereas at any site containing d > 2 particles remain active, and at rate 1, emits two 
particles using independent symmetric random walk steps, leaving (d — 2) particles at the site (observe the 
contrast to the deterministic Abelian Sandpile model, where two particles are emitted in two neighbouring 
directions). Compared to the deterministic sandpile model (see e.g. [7]), much less is rigorously known 
about SSM. It was only recently proved in [26] that there exists pc £ [|, 1] such that the system fixates 
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for fj, < fic and remains active for /x > /Xc- Numerical simulations suggest that Hc ~ 0.9489, it remains a 
major mathematical challenge to prove /Xc < 1. It is reasonable to expect that ARW is a sufficiently good 
(and possibly more mathematically tractable) approximation to SSM and thus captures some of its crucial 
aspects. In particular as A —)■ oo, ARW corresponds exactly to the model studied in [14]. 

ARW can also be viewed as a special case of driven-diffusive epidemic processes. In this process a 
healthy particle does a simple random walk at rate Da > 0 and each infected particle does a simple random 
walk with rate Db > ^, infecting all particles that it steps on, and recuperating at rate A. ARW corresponds 
to the special case Da = 0. This model, in the case Da = Db was introduced by Spitzer in late 1970s, but 
was rigorously studied in detail much later in [15-18]. Although numerical studies have predicted different 
regimes of critical behaviour for Da < Db and Da > Db, not much is rigorously known for the case 
Da ^ Db (see [5] for the study of a related model), in particular, it is not understood whether the behaviour 
for Da = 0 (i.e. ARW) should be the same as the behaviour for 0 < Da < Db- 

For more background and a fuller history of ARW and related processes, the interested reader is referred 
to [12,26] and the references therein; see [24] for a detailed account of non-equilibrium phase transitions in 
lattice models. 

There have been a flurry of rigorous results on ARW in the last few years, mostly in the wake of the 
breakthrough paper [26] , where, as mentioned above, the existence of an absorbing-state phase transition 
was established on Z for both ARW and SSM. This paper crucially uses a construction of the ARW process 
using the Diaconis-Fulton representation and the Abelian property (see the next section). An upper bound 
of one on the critical density for ARW was also established in [26] for any A > 0. The same upper bound 
was also established in [1] for ARW on unimodular graphs using mass transport principle and in [30] for 
general bounded degree graphs using a comparison with certain internal aggregation models. The results 
of [30] also imply that critical density is positive for ARW for any Z'^ for X = oo. This result has recently 
been established for any A > 0 and any d in [31] using a multi-scale argument. As far as the critical value 
is concerned the only result is by [6] for the case A = oo, where it is shown that the system on Z*^ fixates 
for all /X < 1, thus in conjunction with the results of [1,30] establishing /Xc = 1. The behaviour of ARW at 
criticality is also largely open, essentially the only result is for A = oo, in which case the system does not 
fixate at criticality [6]. For more detailed predictions on the behaviour of ARW, see [12]. 

So far we have concerned ourselves with the model where the active particles perform a simple symmetric 
random walk, and that scenario will be the focus of the present work. It is also interesting to consider a 
situation where the random walk steps are biased, and in certain cases, such models are better understood. 
The study of the biased ARW started with an unpublished argument of Hoffman and Sidoravicius (see [6, 
Theorem 1]) that considers the case of totally asymmetric walks and establishes that the critical density 
is /Xc = and further that the system does not fixate at criticality. More recently ARW on Z'^ with 
asymmetric (but not necessarily totally asymmetric) jump distribution has been studied in [28, 32]. It is 
shown in [32] that in d = 1, when the jump distribution is biased. Conjecture 2 holds and Conjecture 1 
holds for small sleep rate. In a very recent paper [28] the same has been established for d > 2 sharpening 
a previous result of [32] which stated that the system does not fixate for small A and /x sufficiently close 
to one. We refer the reader to the lecture notes of Leonardo Rolla [27] for more details on the previously 
known results on ARW. 

Despite this impressive progress in the study of biased ARW, Conjecture 1 and Conjecture 2 have so far 
remained open for symmetric ARW (we shall always refer to the symmetric case by ARW unless otherwise 
mentioned), and it appears that the methods employed in these works cannot be adapted to understand Con¬ 
jecture 1 and Conjecture 2 for symmetric ARW on Z. This is the main contribution of the present paper. 

The proof considers truncated versions of ARW on a finite universe. This technique has been used in 
several other arguments regarding this model in the literature . The new ingredients include a novel use of 
the Abelian property (see Section 3) which allows us to study a slightly different labeled variant of ARW 
which is stochastically dominated by the actual ARW. In other words we construct a coupling of the two 
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systems where fixation of the latter system implies fixation for the former. More details appear in Section 

2 . 1 . 


2. Formal definitions and setup 

Recall that ARW is the following continuous time interacting particle system on Z. Start with particles 
at every site on the line by sampling i.i.d. from any distribution with mean Each particle can be in one of 
the two following states A (active) and S (sleepy). Initially all the particles are in active state. Each active 
particle does a continuous time nearest neighbour symmetric random walk on Z at rate 1. Sleepy particles 
do not move. Also each active particle undergoes the transition A —)■ S at rate A > 0 independent of 
everything else. Sleepy particles undergo the transition 5 + A —)• 2A instantaneously, i.e., a sleepy particle 
at X G Z becomes active instantaneously when an active particle visits the site x. Also the transition A —)• 5 
is observed only if at that time the active particle was the only particle at its site, i.e. the instantaneous 
transitions 2A —)■ A + S' —)• 2A is not observed. 

We follow [26] in formally describing the set up of ARW. Eet p be a formal symbol (whose meaning will 
be clear soon). The state space (2 = (N U {0} U p)^ is the space of all bi-infinite sequence with elements 
of N U {0} U p. Eor any time t > 0, rjt £ Q will denote the state of the system i.e. r]t{x) denotes the 
number of particles at x G Z at time t. rit{x) = p denotes that the only particle at x at time t is in state 
S (is asleep). Eollowing notation from [26] we formally let |p| = 1 so that irrespective of the state of the 
particles, |pt(x)| denotes the number of particles at site x at time t. Eor notational convenience we define the 
following addition and multiplication operations on N U {0} U p to describe the A + 5 — )■ 2A and A — )■ 5 
transitions. 


p+0=0+p=p 

p + n = n + p = n + l. 
p.l = p 

p.n = n for n > 1. 


Eet for X G Z and p G (2, A(p(x)) denote the number of active particles at site x in configuration p. With the 
above notation formally the process evolves as follows: for each site x, we have the transitions p —)■ Tx^yP 
at rate A(pf(x)) 2 l|j/-a;|=i; and p — Tx^pP at rate XA{pt{x)), where 


{ p{z) + 1 z = y 

p{z)-l z = x (2.1) 

p{z) otherwise. 

Similarly 

f \r \ - f P-Vix) z = X 
Tx,p{p){z) I Otherwise. 

Eet denote the law of the process started from an initial configuration distributed according to v. 
Throughout the rest of the article we will focus on u being a product measure with identical coordinate 
projections which in particular implies ergodicity. 


Einally we need to show that such a process is well-defined even starting with infinitely many particles. 
This can be done using the general theory of interacting particle systems developed in [21]. It can also be 
done using an argument of Andjel [2] which shows that (under some mild finiteness condition on v) the 
process is well-defined and can be approximated in a suitable sense by its finite truncations. See [26] for 
more details. 

We say that the system has particle density p if E,y(po(0)) = p. Thus as mentioned before ARW(p, A) 
will denote the activated random walk process with initial density p and sleep rate A (note that we choose to 
suppress the dependence of the distribution v and just keep track of the particle density). 

4 


Definition 4. ARW started from any configuration is said to locally fixate if for every x G Z, r]t{x) is 
eventually constant. Otherwise we say that the system stays active. 

2.1. Key ideas and outline of the proofs. To prove Theorem 1 we first approximate the infinite system 
by considering truncated ARW on large but finite boxes (see Lemma 5). The goal is to then show using 
this approximation that, the number of times a particle gets emitted from the origin can be made larger than 
any finite constant with probability bounded away from 0. More specifically we run several rounds of the 
truncated process with growing intervals, with the origin as one of the end points. The intervals are chosen 
to be growing exponentially in size and it is argued that the total number of times particles gets emitted 
from the origin in the round dominates a Ber( j) random variable. Also by construction we ensure that 
the rounds are independent. Thus the total number of particles emitted at the origin up to the round 
dominates a Bin(^, |) random variable. Using Lemma 5 one then argues that since the activity at the origin 
in the finite volume is arbitrarily large the origin does not fixate in the infinite volume system almost surely. 
Thus Theorem 1 follows. This argument is made precise in Lemma 15 in Section 3.3. 

Proving the claim that the total number of times particles get emitted from the origin in the round 
dominates a Ber(|) random variable is done as follows: we show that for a large enough interval with 
not too few particles it is exponentially unlikely (in the size of the interval) for the finite volume system to 
stabilize inside the interval without any particle touching the boundary (recall that the origin was chosen 
to be one of the end points of the interval). Thus with probability roughly 1/2, particles hit the origin 
(probability of hitting either of the boundary points is close to 1, and both are equally likely to happen by 
the underlying symmetry of the process). 

The technical core of this paper consists of proving the exponential upper bound on the probability of the 
event mentioned above. Roughly we show the following which is made precise as the content of Lemma 
16. Consider a large interval [—r, r], with at least /rr particles (recall that /i is the particle density). We 
show that for A small enough it is exponentially unlikely (in r) that the process stabilizes within the in¬ 
terval [—r, r] (without any particle hitting either boundary). The proof includes a novel use of Abelian 
property which allows the reduction to a study of a labeled variant of the ARW making things tech¬ 
nically convenient (see Section 4). We argue by re-normalizing space: i.e. consider the lattice points 
iTZ = {..., —2K, —K, 0, K, 2K,...}, for AT » ^ and then show at the end of the stabilization process 
for A small enough it is unlikely to have more than a particle asleep in any interval of length K. Since 
we start with enough particles this implies that some of them must escape through the boundary with large 
probability. 

A statistic of fundamental importance in our analysis is the ‘odometer function’ : the number of times 
a particle was emitted from a site until stabilization, (see (3.5) for a formal definition). For other uses of 
the odometer in the study of Abelian systems, see [20]. Using the odometer the remainder of the proof 
proceeds as following: fix a sequence z = {0 = Z-r, ..., Z- 2 k, Z-k, zq, zk, Z 2 K, ■ ■ ■ Zr = 0}. We bound 
fhe probability that the odometer function at the points iFZ n [— r, r] is the sequence z after the truncated 
process on [—r, r] has run until stabilization (this happens in finite time almost surely due to the finiteness 
of the universe). Since by choice Zr = Z-r = 0, on the the above event, —r, r do not emit particles. The 
probability bound as a function of z turns out to be small enough to allow us to take an union bound over all 
such sequences z and end up with e~^'^. 

The above argument relies heavily on a new and novel renormalized variant of the standard Diaconis- 
Fulton representation of Abelian systems (see [8]). One way to run the ARW dynamics, as done in [26], is to 
start with a probability space where at each site on Z one has ‘stacks’ of left, right or sleep instructions used 
by particles when they are emitted, (see (3.2)). The odometer function at a site is then the number of elements 
from the corresponding stack used until stabilization. For our purposes we introduce ‘renormalized’ variants 
of the above where the stacks are only at the renormalized lattice points FfZ and instead of being steps of 
length one, are lazy random walk paths stopped on hitting the nearest point in the set ATZ. The laziness of 
the walk corresponds to the sleep instructions. See (4.3) for precise definitions. 
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2.2. Organisation of the paper. In Section 3 we collect all the basic preliminaries about the truncated 
(finite universe) ARW and the Abelian property. In Section 3.3, we show how to complete the proof of 
Theorem 1 assuming the statement of the key Lemma 16. In Section 4 we develop a labeled variant of the 
ARW dynamics and reduce Lemma 16 to Lemma 18 about the Labeled ARW dynamics. Section 5 develops 
a few technical preliminaries needed for the proof of Lemma 18, which is completed in Section 6 using a 
careful union bound over the possible odometer values. We finish with some open questions in Section 7. 

3. Finite volume dynamics and Abelian property 

For our purposes we now define a truncated version of the ARW(/r, A) dynamics. Recall that the initial 
configuration ryo is distributed according to some product measure v with particle density /r. Let vm be the 
restriction of v on the finite box [—M, M] and := be the corresponding law of the process. The 
next result shows that to prove Theorem 1 it suffices to consider the measures P^. 

To see this consider the following modification of the ARW dynamics. Given a configuration ry on Z, 
and M G N, let rjM denote the restriction of rj on [—M, M]. With the initial configuration rjM, consider 
ARW process with the restriction that particles that move out of [—M, M] are deleted from the system 
(alternatively they fall asleep immediately irrespective of everything else). Since initially there are (almost 
surely) finitely many particles in the system, and since any particle not falling asleep will eventually exit 
[—M, M] it is easy to argue that this process will stabilize almost surely, i.e., all particles will fall asleep in 
finite time. 

The object of interest at this point is the number of transitions happening at the origin, we call this um{0) 
(formal definitions appear in the next subsection). One now considers a coupling all the truncated ARWs 
on one common probability space. We denote the underlying measure by (formal definition appears 
in § 3.1). As an important consequence of the so Abelian property is that under the above coupling the 
sequence it follows that um{0) is nondecreasing in M. Let 

n(0) := lim um{0)- (3.1) 

M^OO 

We then have the following lemma. 

Lemma 5. [26, Lemma 4] Let be a translation invariant ergodic measure with finite particle density 
Ey(ry(0)) < oo. Then 

F^{the system locally fixates) = < oo) G {0,1}. 

We refer to [25], which is an older arXiv version of [26], for a proof of the above Lemma. As already 
mentioned before the 0 — 1 law is a direct consequence of ergodicity. 

In the next subsection we formally define the coupling of the truncated processes mentioned above and 
precisely state the monotonicity property used to define u{0). 

3.1. Coupling of truncated ARW’s. To construct the coupling of truncated ARWs we shall take resort to a 
Diaconis-Fulton representation [8,13] of this process where the process is implemented through a sequence 
of instructions attached to the sites. The advantage of this representation is the Abelian property, which 
allows one to disregard the order in which different steps were performed in certain settings. The relevant 
elements of the Diaconis-Fulton representation in our context was developed in [26, Section 3], and we 
closely follow their treatment for the rest of this and the next subsection. We start by introducing a series 
of notations. Recall the transitions Tx,y and Tx,p from (2.1) and (2.2). Now consider the following array of 
random variables: 

• • • ^(- 2 , 1 ) ?(-!,!) ^( 0 , 1 ) ^( 1 , 1 ) ^( 2 , 1 ) • • ■ 

= • • • ^(-2,2) ?(-l,2) ^(0,2) ^(1,2) ^(2,2) • • • (3.2) 
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where independent for any x £ Z and j G N and moreover, 




Tx,x-i with probability 
Tx,x+i with probability 




2(A+1) 

2(A+1) 

with probability . 


(3.3) 


We would call the ^(x,j) s instructions at the site x and the underlying product measure 0^. Using these 
instructions one can define a discrete time version of the ARW process in the following way: We start by 
defining the notion of stability of a site x with respect to a configuration rj £ Q. 


Definition 6. For any t] £ Q v/e say a site x G Z is stable if r]{x) = 0 or p, and otherwise we call it 

unstable. 


Thus given a configuration 77 G U at each discrete time step t, one can choose an unstable site x and use the 
first unused element from the stack ,^( 3 ,.) and use it to perform the transition to a configuration p' at time step 
(f + 1). We call such an operation “toppling” at site x. Formally we keep track of the number of topplings 
at every site as a function of time t. Let 

h := {h{x) : X £ Z) (3.4) 

where h{x) denote the number of topplings at x, which we will call the odometer function at x. ht{-) will 
be used to denote the odometer function at time t. For later purposes it will be convenient to keep track of 
both the odometer function ht and the configuration pf To this end define the toppling operation at x acting 
on the pair {p, h) by 

^x{ri, h) = (?(x,fc(x)+i)(^), h + 5x), 

i.e. we use C[x,hj:+i)^ the first unused element of the stack at x to topple, and h increases by 1 at x. 

Definition 7. We say that 4 > 3 ; is legal if x is unstable i.e. p{x) > 1. For any sequence a = (xi, X 2 ,..., x^) 
we define the sequence of topplings at xi, followed by X 2 and so on through until x^ by i.e. = 
^x^, ■ ■ ■ ‘hxi. We now say that a is a legal sequence if ^xt is legal for 4>j_i... <l>i for alH = 1,... A:. 

We abuse notation a little to denote by the odometer function after performing the sequence of toppling 
given by a, i.e. for any x G Z, 

k 

ha{x) = '^l{ai = x). (3.5) 

i=l 

We also define the natural ordering on U and the space of odometer functions: for ry, r) G U we say that 
p > pii p{x) > p{x) for all x G Z. Similarly for odometer functions h, h, we say h > hif h{x) > h{x) 
for all X G Z. We also write (ry, h) > {p, h)ifp>p and h = h. 

The most basic property of the above process is the Abelian property which says that given two sequence 
of legal topplings which result in the same odometer function (see (3.5)), the final configuration is the same 
in both the cases i.e. the order in which topplings are performed does not matter. 

Lemma 8. (Abelian property) Given any two legal sequence of topplings a and ol such that ha = ha', then 

T>a(ry) = 4>«/(ry). 

To see why the above is true notice that for each site x,y £ Z with x / yy any pair of legal topplings, 
one at x and the other at y commutes. The sequence of topplings at a single site x do not commute however 
since the compositions of Tx,p and Tx,x+i clearly depends on the order of composition. However given the 
stacks {C(a;,y)} (see (3.2)) and the function /iq(-) one knows in which order the topplings at a single site x 
occurs. 

Before providing a formal proof of the above we now list some basic facts about the toppling operation 
in the following observation. 

Observation 9. The toppling operation as described above has the following properties: 
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i. If a is a legal sequence for r], then 0 ,( 77 ) depends on a only through ha{-)- This is immediate from 
Lemma 8 . 

ii. ‘ho ( 77 ) is non-increasing in ha { x ) and non-decreasing in ha { z ), z ^ x . This is easy by observing 
that a toppling at x can only decrease the number of particles at x , whereas a toppling at another site 
can only increase it. 

hi. If X is unstable in 77 and 770(2:) > 77(2:), then x is unstable in 770. 
iv. Moreover if 770 > 77, then ‘ho( 77 o) > ^^ 0 ( 7 ?). 

Items hi. and iv. above are obvious. 


Proof of Lemma 8. Let \a\ denote the length of the sequence a. Since ha = ha' clearly |a| = \a’\. 
The proof now follows by induction on |a|. Clearly if |q;| = 1, there is nothing to prove. Otherwise let 
a = (xi, X 2 ,..., Xfc) and a' = ( 771 , 7 / 2 ,... t/a:)- Now if xi = 7/1 we are done by induction, considering the 
sequences (x 2 , • • • Xk) and ( 7 / 2 , • • • Uk) with initial configuration ( 77 ) = ( 77 ). If xi 7 ^ 7/1 notice that 

since/icj'( 7 / 1 ) > 1 and/ iq, = there exists a j < fesuchthatXj = 7 / 1 . Let jo be the minimum such j. Now 
consider the sequence = (x^q, 2 : 1 , 2:2 • • • Xj^-i, Xj^+i,..., x^). We now notice that $ 0 ( 77 ) = ( 77 ). 

To see this first observe that is a legal sequence. This is because initially the toppling is legal 
since a' is legal. Also since Xj^ does not occur in (xi,... Xj^-i) by the discussion before the proof, 
commutes with $ 3 ;. 1 < i < jo — 1- Thus ‘^ 0 ( 77 ) = ( 77 ). Now considering aj^ and a' we are done by 

the previous case since both the sequences have the same first element. Thus 

T>«( 77 ) = 03^(77) = 

□ 


3.2. Consequences of the Abelian property. As mentioned before, for our purposes we will be often 
interested in finite ARW dynamics restricted to an interval. We start with the following definition. 

Definition 10. Let 1/ be a finite subset of Z. A configuration 77 is said to be stable in V if all the sites x ^ V 
are stable. We say that a sequence of topplings a is contained in V if all its elements are in V, and we say 
that a stabilizes 77 in y if every x G 1^ is stable in ^a iv) ■ 

Lemma 11. (Least Action Principle, [26, Lemma 1]) Given a set V, if a, fi are two legal sequences of 
topplings such that /3 is contained in V and a stabilizes rjinV then hp < ha, i-C. all the topplings in /? are 
also needed in a. 

Proof The proof is similar to the proof of Lemma 8 and uses induction on |/3|. Let a = (xi, X 2 , • • • Xk) and 
/3 = iui, y 2 , ■■ - yj)- Since /3 is legal, 7/1 is unstable and hence must occur in a. Let jo be the first index such 
that Xjg = 7 / 1 . Thus = (xjp, xi, X 2 ... Xj(,_i, Xjp+i,..., Xk) is also legal and has the same end configu¬ 
ration as a. The proof now follows by induction by looking at the sequences (xi,... xj^-i, Xj^+i,..., Xk) 
and ( 7 / 2 , • yj) and the starting configuration ( 77 ). □ 

As a simple corollary we obtain the following. 

Lemma 12. ( [26, Lemma 2]) If a and /3 are both legal toppling sequences for 77 that are contained in V and 
stabilize 77 in V, then ha = hp . In particular, $ 0 ( 77 ) = ‘h^( 77 ). 

Proof By Lemma 11, / 1/3 < /iq as well as ha < h 13 and hence we are done. □ 

We now state a simple but crucial monotonicity result that allows us define u(0) in (3.1). First we define 
formally the odometer function appearing in (3.1). 

Definition 13. (Truncated odometer) Given 77 G a set of vertices V, and any legal stabilizing sequence a 
as in Lemma 12 define the truncated odometer, 

Uvf) ■■= Urjyf) := haf). 
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(3.6) 


Note that by Lemma 12, ny(-) is a well defined quantity. We also define fhe final configuralion 

Vocy ■= (3.7) 

which is also well defined by fhe above lemma. For our purposes we will only consider finite boxes i.e. 
V = [—M, M], for some positive integer M. For brevify we define 

um-=U[_m,m]- (3-8) 

Lemma 14. (Monofonicify, [26, Lemma 3]) For verfex sefs V C V' and configurafions rj < rj', 

Uri,V < Uri'y. 

Proof. Lef a and /3 be legal sfabilizing sequences for fhe fwo sysfems. By iii. and iv. in Observation 9 we 
see fhaf a is also a legal sequence for 77 '. Thus we are done by Lemma 11. □ 

Using Lemma 14 one can define um- Note also fhaf Lemma 14 implies fhaf does 

nol depend on fhe increasing sefs [—M, M], and any sequence of sefs increasing fo Z would have fhe same 
limif. When fhe underlying configuration is clear from confexf we will denote Uq^oo{-) jusf by u(-) as slated 
in (3.1). A configuration r] is said fo be sfabilizable if u{x) < 00 for all x G Z. (One can easily see fhaf 
almosl surely u{-) is finile everywhere or infinife everywhere). 

We now finish off fhe section wifh fhe slafemenl of some fechnical lemmas and deduce fhe proof of 
Theorem 1 from fhem. The resl of fhis article will be devofed fo fhe proof of fhe lemmas. 

3.3. Proof of Theorem 1. Before proceeding furlher wifh fechnical argumenfs, we now discuss fhe key 
sfeps in fhe proof of Theorem 1. By Lemma 5 if suffices fo show fhe following: for any /i > 0 for small 
enough sleep rate, almosl surely, fhe sequence {riM(0)}MeN (see (3.8)) fakes arbifrarily large values. To 
show Ibis we look af fhe following sequence of inlervals of exponenlially growing size. 

= [ 0 , 2 ^] 

^2 = [ 0 , 28 ] 


U, = [ 0 , 23^+2] 

Recall T] is fhe initial configuralion on Z disfribuled as fhe producf measure v. Lef denofe fhe reslricfion 
of rj fo fhe interval li = [2^^, 3 • 2^^]. Observe fhaf Ifs are jusf fhe middle halves of fhe inferval and 
are mulually disjoin! for £ € Z. Since 1 / is a producf measure, are independenl samples from 

measures respectively where denotes fhe coordinale projecfion of fhe producf measure u onlo 

fhe interval 

We Ihen consider running rounds of Iruncaled ARWs wifh initial dala for and verfex sel Vg for 
i = 1,2,.... Thai is, in fhe £^^ round, we ignore all fhe particles fhaf are nol initially in and slabilize fhe 
initial configuration on fhe inferval Vi using fhe elemenfs of fhe slacks (see (3.2)) fhaf have nof been 
used upfo {£ — 1)*^ round. Notice fhaf in fhis sequential process, we are ignoring cerlain active particles 
from previous rounds, however using Lemma 14 fhis will suffice for our purpose of oblaining a lower bound 
on fhe odomeler counfs. Lef we denofe fhe increase in fhe odometer counl al 0 during fhe £^^ round of fhe 
process described above. We shall show fhaf ff^i we becomes arbifrarily large. To fhis end we have fhe 
following lemma. 

Lemma 15. are mulually independenl. Furlher, for /i > 0 and A sufficienlly small depending on 

/i we have fhaf for all large £, 

u (8) £^{wi > 0) > -. 

We posfpone fhe proof of Lemma 15 for fhe momenl and inslead show firsl how fhis implies Theorem 1. 
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Figure 1. The interval of activity Vi in the round. 

Proof of Theorem 1. It follows from Lemma 15 that = oo, u (S> ^-almost surely. Observe that at 

the end of the round, the odometer function at the origin takes value Now recall uv^(O) from 

(3.6). Thus by Lemma 5 we would be done once we show for £ > 1, 

i 

uvt{0)>^Wi. (3.9) 

^=1 

To see the above notice the following: all the steps up to the £^^ round are legal even though they need not 
be stabilizing in Vf. Indeed while running the rounds we ignore particles from previous rounds which can 
affect the process in two ways. In our sequential operation described above, a particle that is asleep at the 
end of a round is static throughout all the next rounds, whereas while stabilizing rj on these particles can 
be woken up by movements of particles in later rounds. Also some of the sleep instructions implemented 
during the £*^ round might not be implemented in the actual stabilization process on Vi. This is because 
other particles from previous rounds might be present at the site (recall that these particles are ignored while 
running the £^^ round, so in the round these sleep instructions might actually have been implemented). 
However crucially, even though the sequence of topplings up to round £ need not be stabilizing, all the steps 
are legal with respect to the stabilization process on the set Vi with initial particle configuration being r] 
restricted to V^. Thus (3.9) follows by Lemma 11. □ 

It remains to prove Lemma 15. The main step is to show that if the number of particles in the configuration 
is not too small, it is exponentially unlikely in the length of the interval Vi that during the £^^ round of 
the operation described above (i.e., stabilizing the initial configuration of particles in the interval Vf) 
none of the particles actually move out V^. Formally we define 

Ar ■■= {max{M,^j_2,.,2r](-2r'),M*J_2r,2r](2r)} = O} (3.10) 

In fhe following lemma we bound fhe probabilify of Ar- Its proof will take the rest of the paper. 

Lemma 16. Fix /x > 0. Let be a particle distribution supported on [—r, r] such that the total number of 
active particles in tjr is at least /xr. Let denote the law of the ARW dynamics with initial configuration 
and the stacks with law (as in (3.2)). Then for A sufficiently small depending on /x then there exists 
c = c(/x, A) > 0, such that for all r 

(A) < 

where u\^- 2 r, 2 r] (') denotes the odometer function (see (3.8)). 

We finish off fhis section by proving Lemma 15 using Lemma 16. 

Proof of Lemma 15. Recall fhaf since fhe infervals are disjoin!, fhe inifial configurafions are indepen- 
denf samples from fhe disfribufions for £ = 1,2,.... Also observe fhaf, because of fhe independence of 
the elements of the stack, the joint distribution of the unused elements of the stack at the end of the {£ — lY^ 
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round (i.e. the stack that is being used for the round) is the same as the law 0 ^ of the original stacks, 
and further this is independent of the elements of the stack that have been used in the first {£ — 1) rounds. 
This implies that are mutually independent and further that for each £>1, the distribution of W(, (under 
V ® £ 0 ’) is the same as the law of uy^ (0) under vi^ ® 0 ^. 

Let denote the number of particles in the configuration distributed as v^,. Due to an obvious trans¬ 
lation invariance in the system, Lemma 16 implies that 

® 00 (^max{riy«(0),uvj(2^^+^)} = 0, 

Observe that by strong law of large numbers for £ sufficiently large we have > /u2^^) > 0.9 and 

hence for £ sufficiently large 

® ^ (^max{ttyj(0),Myf(2^^+^)} > 0^ > 0.9 - 

Notice that the interval le and are both symmetric about 2^^+^. Because is a product measure with 
identical marginals, and the dynamics of ARW possesses an inherent left-right symmetry, it follows by a 
reflection about 2^^+^ that uv^{£)) and uv^(2^^+^) has identical distributions under the measure ® £ 0 *. 
The lemma follows. □ 

The rest of this paper is devoted to the proof of Lemma 16. Recall the stacks and the underlying measure 
00 from (3.2). In case of finite vertex sets (which will be the focus of our analysis) by Lemma 8, the 
final odometer function and the particle distribution is just a function of the stacks and the initial particle 
configuration. However for our purposes we need a different set of stacks, labeled particles, and a certain 
toppling rule. This is described in the next section and we refer to this as Labeled ARW dynamics. Thus to 
distinguish the two, we denote the description of ARW so far (equivalently the measure space corresponding 
to J' (see (3.2) and the initial particle configuration rf) as Unlabeled ARW dynamics* The next section also 
provides a coupling between the two ‘processes’. 

4. Labeled ARW Dynamics 

We start by outlining how we prove Lemma 16. Remember that the lemma states that if we start with a 
typical particle distribution in [—r, r] (and 0 outside this interval) then it is unlikely that when we stabilize 
the odometer at both —2r and 2r is 0. Typically we shall have about 2|Ur many particles in our system 
and with very high probability we will have at least /rr particles in the interval [—r, r]. Notice that on the 
event Ar, which was defined in (3.10), all the particles eventually fall asleep in the interval [—2r, 2r]. Now 
consider K a large integer Divide the interval [—2r, 2r] into subintervals [iK, (i + 0)K] of length 

K each. In the final configuration, we typically expect only a small number of particles asleep on each 
interval [iK, {i + 1)K] if A is sufficiently small depending on K. Indeed, whenever a particle is released 
from the site iK it is extremely likely (by taking A small) that this particle will reach either (i — 1)K or 
(i -|- 1)K without falling asleep in between. Further, in such a case, it wakes up all the sleepy particles in 
either [iK, {i + 0)K] or [(i — 1)K, iK], As A is large it is highly likely that all of these recently woken 
particles will make it to either {i — 1)K ox {i + 1)K before falling asleep again. Hence it is unlikely that 
in the 4r/iT many intervals of length K on average there will be iT/r/4 particles asleep (since K 3> /r“*). 
Thus not all particles (at least fir) will fall asleep before at least one particle makes it to —2r or 2r. 

To make this intuition into a proof, we shall introduce a different set of stacks containing instructions, 
where, instead of a single step, the instructions now will consist of a random walk path that tells the particle 
(started at iK) its steps until it reaches iK ± K. 


*Note the abuse of notation since there is no canonical process in the discrete time setting and any order of toppling eventually lead 
to the same statistics by Lemma 8. 
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4.1. Random Walk Stacks. Let A > 0 and consider the lazy symmetric random walk S'a(-) with laziness 
on Z started from the origin. By this we mean 5^(0) = 0 and in each discrete time step, f = 0,1, 2,... 

S^{t + 1)-Sx{t) = < 


independent of everything. Now for any positive integer K let t{{—K, K}) denote the hitting time of the 
set {—K,K} for the walk S'a(-) and let 

Sk,x{') = 'S'A(min(-,r({-it:, K}))) 

be the killed random walk, i.e. Sx{-) stopped on hitting {—K, K}. 

We now define a probability space where instructions will be the killed random walk paths instead of 
single steps as in (3.2). We consider the renormalized lattice, 

KZ = ...-2K, -K, 0, K, -2K, (4.2) 

and the following set of instructions at the points iLZ (which we shall sometimes refer to as lattice points): 

• • • C(-2,l) C(-l,l) C(0,1) C(l,l) C(2,l) • • • 

= • • • C(-2,2) C(-l,2) C(0,2) C(l,2) C(2,2) • • • (4.3) 

where C{x,j) are independent for any x G Z and j G N and moreover, 

C{xj) = xK + Sk,\{-), (4.4) 

i.e. the elements C{x,j) are i.i.d. copies of random walk paths Sk,x{-) started at xK and stopped on hitting 

xK ± K. We will still continue to think of the random walk steps to be distributed as the elements in 

(3.2) namely: 

Tz,z-i with probability 

T^^z+i with probability (4.5) 

Tz,p with probability . 

where z is the current location of the random walk. That is left and right steps in the random walks would 
be interpreted as instruction to go left or right, whereas lazy steps would be interpreted as sleep instructions. 
However, how these instructions act upon the configurations will be slightly different from the standard 
ARW dynamics of § 1 as explained below. We call the above product measure 3^^ analogous to 3^. 

4.2. Labeled Particles and Labeled ARW Dynamics. With the above stacks we now define a variant of 
the ARW dynamics from the discussion following Definition 6. To define fhe modified dynamics, we shall 
resfricf ourselves fo initial particle configurations supported on iTZ. In this setting the particles will be 
labeled. The labels will keep changing over time and will be from the set Z x N. 

Labeling Scheme: Suppose we start with a configuration rj^ supported on ATZ. The first co-ordinate of 
the label of all the particles starting at iK will be i. The second label is arbitrarily chosen to be: 1 to the 
number of points on that lattice site. That is, if there are n* particles at iK to start with they are labeled 
(i, 1), (i, 2),..., (i, rii) in an arbitrary way. See Figure 2. 

As mentioned above the labels are not static but evolve with time. Any particle of label {i,j) will change 
its label when it hits the lattice points {i — l)K or (i + 1)K in which case the label would change to (^^ j) 
where i' = i 3 I according to the lattice point hit and j is the smallest positive integer such that the label 
{i',j) has never been used before in the history of the process. A particle with label {i, •) will sometimes be 
referred to as a particle emitted from iK. 

We now introduce a toppling scheme for the labeled process analogous to § 3.1. To do that we first 
associate with each particle a random walk from the stacks (4.3). 
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Figure 2. Initial labeling 


Associating Random Walks with Particles: The first time a particle acquires label (a, b) (recall that our 
labeling rule implies that at this time the particle must be at aK), it gets associated with the random walk 
C{a,b) which stays associated with the particle until it changes its label next. 

Next we define how to use the random walks to topple particles. Informally we do the following. When¬ 
ever we decide to topple a particle with label (a, 6), we use the first unused step of the random walk C{a,b)- 
Let r be the total number of steps in C{a,b) until the random walk hits aK ± K). Whenever the particle 
is toppled later before r, say for the time, the toppling step is the step of C(a,6)- We now provide a 
formal description. 

Toppling Rule and Labeled ARW dynamics: Recall the transitions A —)• 5 and A + S' —)■ 2A between 
particles in active and sleepy states from § 3.1 and that only A particles could be toppled. We introduce a 
particle interaction which is a slight variant of that dynamics, where the particles with different first labels do 
not interact. We shall call this variant, the Labeled ARW dynamics and it will always be defined on a finite 
universe V. That is, we shall fix a finite subset V of Z and and particles outside V become instantaneously 
inactive and stay inactive for ever. To distinguish from the standard ARW dynamics from § 3.1, we denote 
the state of active and sleepy particles in the Labeled ARW dynamics by A and S respectively. 

Consider stacks distributed according to and a particle configuration r/* distributed supported on AZ. 
Fix V and let all particles in V initially be in A state. Consider the labeling scheme described above. For our 
purposes we shall work with particular order of topplings. To this end we start by defining the lexicographic 
order on the labels: i.e. (a, b) < {c, d) if a < b or a = b and c < d. 

At each time t, we topple the A particle with the smallest label, say (a, b), which is well defined as we are 
restricted to the finite set V ? To topple, use the first unused instruction of random walk C{a,b)- The left and 
right step instructions from the random walks are interpreted as before and used to move the particle one 
step to left or right respectively. The lazy steps in the random walks are interpreted as sleep instructions. 

However, in the Labeled ARW dynamics the A + S —2 A transition occurs only when the label of the 
two particles have same first co-ordinate. Let us illustrate this with an example. Consider an S particle at 
y G {aK, {a + f)K) with label {x,j). Notice that by definition either x = aorx = a + l. Now if an A 
particle is toppled from ?/ ± 1 and moves to y, the transition A -|- 5 —)• 2A occurs only if the A particle has 
label of the form (x, •). Thus in words the above says that a particle which was emitted from a lattice site 
aK G KTj can be woken up only by another particle emitted from aK. 

Also in the Labeled ARW dynamics the A —)■ S transition is allowed to happen only when the label of 
no other particle at that site has the same first co-ordinate. Notice the difference with A —)■ S' transition. As 
observed before the particles with different first labels do not interact. ^ 


Note that in Unlabeled ARW dynamics, particles were unlabeled and hence the definition of the process involved toppling sites 
rather than particles unlike here. 

As a consequence that now we are allowed to have two S particle per site provided they have different first label. Moreover, if 
there are two S particles at y € {iK, (i + l)K) (one with first label i and another with {i -|- 1)) and an A-particle with first label 
i moves to y. It then wakes up the S-particle at y with label i, but not the one with label i -I- 1, i.e., in such a situation we have a 
A + 2S ^ 2 A + S transition. 
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The system is said to stabilize if all the particles in 1/ are in 5 state (recall that a particle that moves out of 
V is inactive forever). As before it is easy to argue that as V is finite and the number of particles is finite the 
system stabilizes almost surely in finite time. We define the odometer function and the final configuration 
analogous to (3.6) and (3.7) respectively. 

Note in Unlabeled ARW dynamics there was no toppling rule. Hence to ensure that (3.6) and (3.7) are 
well defined we had to invoke the Abelian property (Lemma 8). However in the Labeled ARW dynamics 
the toppling rule is fixed and hence there is no ambiguity in the definitions of the final odometer count and 
the final configuration. We denote the odometer function and the final configuration by 

ri*y(-) and 7/*,oo,t/ (4.6) 

respectively. Note the dependence on the initial configuration is suppressed. 

Now we want to relate the Labeled ARW dynamics with Unlabeled ARW dynamics. Since by definition 
in the labeled process, particles with different first labels do not interact, one would expect more active to 
sleepy transitions and less sleepy to active transitions, than the unlabeled process. So started with the same 
initial configuration it is expected that the Labeled ARW dynamics is quicker to stabilize (in a stochastic 
sense). To make this formal we need to define a coupling between the measures ^ and (the law of the 
stacks used in the unlabeled and the Labeled ARW dynamics respectively, see (3.2), (4.3)). 

4.3. Coupling Labeled ARW Dynamics with Unlabeled ARW dynamics. Fix a finite subset U C Z. 
Given an initial particle configuration ry supported on V, we define a natural coupling of the stacks 
and having laws y and respectively. The coupling is run in three rounds. Note the labeled process 
is only defined when the initial configuration is supported on KZ. This is what the first round achieves. 

First Round: We run Unlabeled ARW dynamics in V where we topple all active particles in Z\KZ. 
After the end of this round the system only has active particles at U n KZ. Call the random configuration 
of active particles fj, supported on KZ. The distribution of ry is a function of ry. 

Second Round: In this round we describe a ‘natural’ coupling of Unlabeled ARW dynamics and the 
Labeled ARW dynamics starting from ry. The coupling V will be Markovian and will ensure that the ARW 
processes are exactly the same for all times. This is done in the following way. Recall the labeling of the 
particles and toppling rule described in Section 4.2. Also recall that by the Abelian property (Lemma 8) 
the final statistics of Unlabeled ARW dynamics are independent of the order in which particles or sites are 
toppled. Thus in order to couple with the labeled process we choose the same toppling rule as the latter. 
Namely, we topple the same particle in both the processes and couple the stacks y and y^ such that the two 
processes stay identical throughout. Note that this is possible due to the following two reasons: 

• The steps of the random walk .) in (4.3) have the same distribution as .) in (3.2), as already 
remarked in (4.5). 

• Arguing inductively, if the particle configuration stays the same in both the versions under V up to 
the toppling, since an A particle is always an A particle, the [t + 1)*^ toppling would be legal in 
both the processes. 

As a consequence of the above the stabilizing sequence in the Labeled ARW dynamics would always be a 
legal sequence in Unlabeled ARW dynamics. Thus this round of the coupling runs until the Labeled ARW 
dynamics has been stabilized in V. 

Note that potentially many of the particles in rj that were not in fj and hence not toppled in the second 
round could have woken up in the process rendering the configuration running in Unlabeled ARW dynamics 
unstable. The last round takes care of these remaining particles. 

Third round: Independently sample the rest of the stacks y and y^,. Use the former and run Unlabeled 
ARW dynamics until stabilization toppling arbitrarily. By Lemma 8 does not change the final distribution. 

Thus we have completed the description between the two processes. The following observation is trivial 
from the definition of P^. 
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Observation 17, Under the coupling described above, uv{-) started from rj is lower bounded by 
(see (4.6)) starting from fj. 

To prove Lemma 16 it suffices to lower bound ri* j_ 2 r, 2 r] (') started from fj. Remember that r) is a random 
particle configuration supported on iTZ whose distribution depends on tj, but is independent of the stack 
and rj is an initial configuration supported on [—r, r] with at least /rr particles. This is the content of the next 
two lemmas. 

Lemma 18. Fix /x > 0. There exists K = K{^) and A = A(/x, K) > 0 such that the following is true for 
all r. Consider an initial particle configuration rj supported on [—r, r] n iTZ with at least ^ particles. Then 
there exists c > 0 such that for all r 

(4.7) 

where =^* 5 ^ denotes the law of the Labeled ARW dynamics started with the initial particle configuration rj 
and Ar was defined in (3.10). 

Lemma 19. Recall the coupling Vr/ from above. For any rj supported on [—r, r] with at least /xr many 
particles, for A sufficiently small there exist c > 0 such that 



where |r)| denote the number of particles in fj. 

Proof. For any particle in t/, at y G ((x — l)iF, xiT), notice that for A sufficiently small depending on K, 
each particle has a chance at least | of reaching {(x — 1)K, xK} before falling asleep and these events are 
independent of each other. Since all such particles are in ij, the conclusion now follows from the fact that 
there are at least particles in 77 and a large deviation estimate. □ 

We finish off this Section with the proof of Lemma 16. The rest of this paper is devoted to proving Lemma 
18. 

Proof of Lemma 16. Given /x > 0 we get K and A from Lemma 18. The proof follows from a combination 
of Lemma 19, Lemma 18 and Observation 17. □ 

In the following section we develop the remaining technical results needed for the proof of Lemma 18. 

5. Technical Preliminaries 

Recall the set up of Lemma 18: fix iF = iF(/x) sufficiently large and r sufficiently large depending 
on K and a particle configuration r) supported on [—r, r] n iTZ that contains at least ^ many particles. 
Consider Labeled ARW dynamics on [—2r, 2r] starting with initial configuration rj. We shall denote the 
law of this process by Let 7)00 denote the final configuration and for xK G (—2r, 2r) let fjx^oo 

denote the final configuration of the particles labeled (x, •). It is trivial to observe that rj^^oo is supported on 
((x - l)iT, (x + l)K). 

Also recall that u^^y_ 2 r, 2 r]i.x) counts the number of times a labeled particle is toppled at a site x up to 
stabilization. However for our purposes we need to determine how many elements (random walk paths) 
from the stacks in (4.3) have been used. For xK G iTZ, define the renormalized odometer as: 

M(x) := max{j : a particle labeled (x,y) was emitted} , (5.1) 

i.e. the number of random walk paths C{x,-) iri (4-3) that were used until stabilization. 

Now we define certain statistics of the process which will be used in the proof of Lemma 18. For each 
xK G [—2r, 2r], we keep track of the net number of particles labeled (x, •) (emitted from xK) hitting 
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(5.2) 


(x — 1)K, xK — 1, xK + 1 and {x + \)K respectively by 

L{x — 1 •(— x) = L{x — 1 •<— x){r]), 

L{^ x) = L{^ x){fi), 

R{x —>) = R{x —^)(r/), 

R{x —)■ X + 1) = R{x —)■ X + l)(r 7 ) 

respectively. Note that we are counting the net number of particles in all the counts, i.e. if a particle moved 
back and forth across xK a few times to eventually land somewhere in the interval {xK, (x + it just 

contributes 1 to R(x —)■) but not to R(x —)■ x + 1). On the other hand he particle since and 0 to others. See 
Figure 3. Next we make a crucial observation that the statistics in (5.2) are “local” functions of the Labeled 



Figure 3. The particle (x, 1) eventually lands to the right of xK and hence contributes one to 
R(x —>). On the other hand the particle (x, 2) hits (x — 1)K and hence contributes one to L(<— x) 
and L{x — 1 t— x). 

ARW dynamics. This is the reason behind defining the Labeled ARW dynamics in such a way that particles 
with different first labels do not interact. We make things formal below. 

5.1. Single site dynamics. For any x G Z consider Labeled ARW dynamics in 14 = ((x — 1)K, (x + l)K) 
started with m particles at x. Then the quantities analogous to (5.2) will be denoted by 

(m), (m), R^^ (m), (m) (5.3) 

respectively. Thus the above can be thought as a localized version of the Labeled ARW dynamics. The final 
configuration in this process will be denoted by 

Consider Labeled ARW dynamics started with initial configuration fj stabilised in (—2r, 2r). The parti¬ 
cles with label (x, ■) are toppled only at certain times say n, r 2 ,.... Thus the quantities in (5.2) are functions 
of the configurations at these times. By definition 

{(x, 1), (x,2), • • • , (x,M(x))} 

are the different labels of such particles. 

Now compare the above with the restricted system where the vertex set is ((x — l)i^, (x + l)K) and the 
initial configuration of M (x) labeled active particles at xK. Notice that owing to the labeling scheme and 
the labeled toppling rule the movement at times ri, r 2 ,... as described above is exactly the same as in the 
single site dynamics. All this can be formalized by an easy induction on M (x), hence we have the following 
observation. 

Observation 20. For a fixed initial configuration r) recall the definition of M(-) given in (5.1). For each 
X G Z such that xK G (—2r, 2r), the quantities in (5.2) are the same as the respective quantities in (5.3) 
with m = M(x), where both dynamics use the same stack (see (4.3)). Moreover r]^^{M{x)) = fjx^oo, 
i.e., the eventual distribution of the particles labeled (x, •) is same in both processes. 
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Notice the restricted process clearly is just a function of M(x) and the stack at x( {Qx, )})- Hence the 
above observation says that after the labeled system has fixated the distribution of the particles labeled (x, •) 
and the other local statistics defined in (5.2) are just functions of the stack {Qxj) : j £ N} at x and M{x) 
(the renormalized odometer value at x). We also note that there might be particles labeled (x ± 1, •) in the 
interval ((x — l)K, (x + l)i^) after fixation. However the stack at x allows us to determine the distribution 
of particles labeled (x, •). This fact will be crucially used in the proof of Lemma 18. 

6. Proof OF Lemma 18 

The proof follows by a careful union bound over the possible values that the function M(-) (see (5.1)). 
We start by describing the rough idea before providing formal arguments. Recall the event from (3.10). 

Ar := {max{t6*j_2r,2r](-2r-),u*j_2r.,2r](2r)} = 0 } ( 6 . 1 ) 

and that the initial configuration fj is supported on [—r, r] with 2/ir > \ fj\ > Note that in Lemma 18 we 
only assume the lower bound on \ fj\. However by monotonicity (Observation 9 iv.) it suffices to assume the 
upper bound is the same as the lower bound. 

Observation 21. On Ar, no particle hits {—2r, 2r} and hence at least ^ particles fall asleep inside 
(—2r, 2r). 

We shall need the following notation to keep track of the net flux of particles across lattice sites. Namely 
for every x G Z, let 

= R{x —^) — L(x X + 1) (6.2) 

F~ = T(-(— x) — R{x — 1 —)> x) 

where R{-) and L(-) are defined in (5.2). Thus F~ and F+ denote the net number of particles moving across 
the points xK — .5 and xK + .5 respectively. 

For a positive integer oq and integers /q , f^, let V = 'D{ao, fQ , /q ) denote the event that the following 
equalities hold. 


M(0) = ao, F+ = f+,F^ = /q-. (6.3) 

For a fixed triple (uq, /(7) we shall show that V n Ar is unlikely, and Lemma 18 shall follow from a 
union bound over different triples. It turns out that we need to consider two cases, we start with the following 
case, where M(0) = oq is ‘small’, i.e. bounded by a polynomial in r. 

6.1. M(0) is small. Recall the probability measure from Lemma 18. 

Lemma 22. Fix {qq, , Jq) as above with oq < r® and |/(j*"|, | /()”| < 2^r. There exists K = K{u) 
sufficiently large such that for A sufficiently small and for some constant c = c(/r, A) > 0, for all large 
enough r, 

n T>(oo, Jq, /o")) < (6.4) 


For the rest of this subsection, fix (oq, /(j*", /q ) as in Lemma 22. Recall also that t), the initial 
configuration of active particles, supported on iTZ n [— r, r] has already been fixed. 


Given the stack at 0, (i.e. {C(o,j) • j ^ ^1) using Observation 20 and the fact that M(0) = oq, one can 
compute ro := i?(0 —)•). Since Fq = by definition (6.2), 

4 := L(0 ^ 1) = ro - f+. 
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(6.5) 


Now li, together with the stack at K, i.e., {Cij : j £ N}, can be used to calculate M(l) upto some 
(typically small) error. In general, given £x, one can work out the value of M(x) upto some error. Let 
£x{-) '■= where the latter was defined in (5.3). We have the following lemma. 

Lemma 23. For any x G Z and every realization of the stack at xK (i.e. {C(x,j) ^ j £ (see (4.3) ) £x{-) 
is a non decreasing function with jump size at most 2K. 

Proof. Notice that for j > 0, we have 

^x{j + 1 ) = ^x{j) + (■* 

where £^ denotes the net number of particles going to (x — 1)K in the following process: Start with ini¬ 
tial configuration of S particles distributed according to r]^^{j). Add an additional A particle with label 
{x,j + 1) at xK (this will wake up S particle present at xK, if any), call this initial configuration r/*. Now 
run Labeled ARW dynamics with r/* stabilizing on ((x — 1)K, (x + l)K) using the remaining parts of 
corresponding stack elements. This immediately shows monotonicity. To see the bound on the jump length 
note that < | 7 y*|. Since 77 * can have at most two particles at xK, and at most one particles at any other 
site in ((x — l)K, (x + 1 )K), the second assertion of the lemma follows. □ 

Now it might turn out that for some realization of the stack {C(i,j) ■ j £ ^i(') does not take the value 

ro — /q and hence 'D{ao, /J*", ff) cannot hold. On the other hand since the function £i{-) is not strictly 
increasing it might also happen that £\{-) = tq — on a non empty interval of integers. Let the interval be 

[M*(l), M*(l) + t(;(l)], ( 6 . 6 ) 

for some positive integer M*(l) and ui(l) > 0, i.e. for all y G [M*(l), M*(l) + t(;(l)], £i{y) = xq — Jq. In 
this case, based on the information so far (i.e. V{ao, Jq , ff) and the stacks at 0 and K) one cannot determine 
the exact value of M(l) and it could potentially be any element in the interval [M*(l), M*(l) + m(l)]. 

Now, observe that if we knew the value of M(l) G [M*(l), M*(l) + r(;(l)], it would also be possible to 
work out the value from the stack at K and hence the value of £2 as well. So recursively, assuming we 
know that value M(x) and using stack at (x + l)K, using Lemma 23 and Observation 20, we can identify 
a interval [M*(x + 1), M*(x + 1) + w{x + 1)] in which M(x + 1) must be contained. Such an interval 
being empty implies our assumption on the values M(0), M(l),..., M(x) (and also F^) must be incorrect. 
Formally we use the following parametrization. 

Definition 24. Let be the number of labeled sleepy particles in the configuration r]x^{j) (by 

definition all of them have label (f, •)) between {Ki,K{i + 1)) and similarly Si^i-i{j) denotes the number 
of sleepy (f,-) labeled particles between {K{i — l),K{i)]. Also following Observation 20 let = 

and Si^i-i = Si^i-i{M{i)) be the (i, •) labeled sleepy particles in fjx^oo- 

Fix a sequence of non negative integers, 

f 2 r 2 r 1 

- ^ ■ • ’?f(“ 2 ),y(-l),y( 0 ) = 0 ,y(l), 2 /( 2 ),.. ■ , 2 /(^ ~ 

Given the triple (oq, f^, ff), we now define a filfrafion {^i}i>o generafed by fhe variables Q..) (also 
depending on fhe initial particle configuration f) and y) as follows. For convenience of reading, we shall 
denofe fhe filfrafion by fhe sef of random variables fhaf generafe if. We sfarf by defining, 

^0 := {C(0,1)) C(0,2)) • • • C(0,ao)}' (6.8) 

We now recursively define a sequence of funclions all of which will be progressively measurable wifh 
respecf fo fhe filfrafion (yef fo be fully defined). For f > 0, assume £ 1 ,... ,£i \ ai,... ,ai and 
have already been defined. Lef 

4+1 := i?f4(a.) - ft 
18 


(6.9) 


We now define the filtration recursively for y such that y{0),y{l ),..., y{i) is realizable (also recur¬ 
sively defined below) . Assume have been defined. Analogous to (6.6), define M*(f + 1) and 

w{i + 1) by 

M*(f -f 1) = min{j > 0 : = ^i+i}; M*(i -f 1) + w{i -(- 1) = max{j > 0 : = ii+i} 

if such numbers exist. Call the sequence {y{0),y{l),.. .y{i + 1)} realizable if (a) either z = 0 or 
{?/(0), y(l),... y{i)} is realizable and (b) M*(z + 1) and w{i + 1) as above exist and y{i -f 1) < w{i + 1). 
In such a situation define 


®i+i — M*{i -l-1) -l- y{i -f 1). 
For a realizable sequence {y(0), y{l ),... y{i + 1)}, define 


( 6 . 10 ) 



Now by definition i?j+i_^(aj+i), 5i+i^i+2(oi+i) are ,^i+i measurable. Also define 

/j^l := '>l{iK+l,{i+l)K] + 


( 6 . 11 ) 


where |f)s| denote the number of particles of fj contained in B C Z. That is ^i+i measurable is now 
a simple consequence of (6.11). Thus by the above discussion, using (6.9),(6.11) and (6.10) repeatedly we 
inductively construct the filtration for realizable sequences. 

Now for z < 0, we construct in the exactly similar way a filtration and a sequence of functions 
w{i), ai, ri and f~, recursively going from right to left and starting with (ao,/g”). The se¬ 
quence {y{i), ■.. ,y(0)} is called realizable if y{i) < w{i). Finally, we call y realizable if {y(—^ + 
1), • ■ •, 2 /(- 2 ), y(-l), y(0)} and {y{0),y{l),y{2 ),y{^ - 1)} are both realizable. 

Following is a consequence of Observation 20: suppose (z/(0),..., y{i)) (resp. (y(z),..., z/(0))) is real¬ 
izable, and M(z) = ai- Then, on D(ao,/()^,/()")n.4r, there exists z/(z-|-l) > 0 (resp. z/(z —1) > 0)suchthat 
(y(0),..., y{i + 1)) (resp. (y(z — 1),..., z/(0))) is also realizable. This leads to the following observation. 

Observation 25. On 'D{ao, , Jq) r\ Ar, there exists a (unique) y as in (6.7) such that y is realizable and 
M(z) = M^{i) + y{i) for each z. Further for such an y, we have ij = L{j — 1 •(— j) and for all 

j > 0 where ij and are defined as in (6.9) and (6.11) respectively. 

It might be instructive to think of the path y as a path in a tree of possibilities of the odometer values at 
various points, out of which the process picks only one. 

Denote by TZ{y) the event that y is realizable and M(z) = M^{i) -|- yi for all z. Also for j > 0 (resp. 
j < 0), let Tl{y;j) denote the event that {y(0),... y{j)} is realizable (resp. {y{j), ■ ■ ■, y(0)} is realizable) 
and M(z) = M^{i) + yi for all z G {0,1,..., j} (resp.ll z G {y, j -|- 1,..., 0}). We now state for any 
y and the triple (ao, fo , f^) as in Lemma 22, a quantitative upper bound on the probability of the event 

h{y) n D n A- 

Lemma 26. There exists c > 0 such that for any K sufficiently large depending on y, there exists \k 
sufficiently small, such that for A < Xk, and y as in (6.7), 


'^*,^(7^(|/) n D n A) < 


To start with recall from Observation 21, on Al, 



( 6 . 12 ) 


Fix y as in (6.7). To prove Lemma 26 we prove the following. 



(6.13) 
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The proof of Lemma 26 then follows by Markov’s inequality and the observation that on the event 
We now proceed toward proving (6.13). We shall prove the following for all j G {1,..., 2r/iT 


E 


g[4j + l+4i-l]l [7^(y; j) r^Vr^A\ I 


Similarly for y G {—2r/iT + 1,..., —1}, 


E 


g[4j + l+4i-l]l [7e,(y; j) n n I + 1 






vr\A 

(6.14) 
- 1 }: 

(6.15) 

(6.16) 


Clearly the above bounds for all j implies (6.13). Note that in the above equations we do not consider 
*50,1 + *50,-1 • However we naively bound it by 2K (since at most 2K particles with first label 0 can fall 
asleep in the interval (—iT, K) ). This accounts for the extra term in (6.13). We shall only show (6.15), 
the other bound will follow similarly. 

Before proceeding further we need some notation: recall the random walk paths q in (4.3) are lazy 
(with parameter A) stopped on hitting the nearest renormalized lattice point (multiple of K) to the left or 
to the right. By symmetry each random walk path has probability 1/2 of hitting either neighbour. Call 
a random walk path q a left instruction (resp. right instruction) if it is killed at hitting (a: — 1)K 
(resp. {x + 1 )K) before taking any lazy step. Clearly for any x > 0 A = X{K,x) sufficiently small 
is ^ instruction) > ^ — x and similarly for a right instruction. Fix any e > 0. Let R be such 
that the probability of having a right instruction in R i.i.d. samples of the random walk paths is at least 1 — e. 

For every x G Z, given we define the following stopping times measurable with respect to the filtration 


^xii) {C(a;,l)) ■ • • C(x,t)}- 

Recall lx{-) from Lemma 23 and define 

Ti = ri(x) := infjf : ^x{t) > (-x — 3iTi(} 


(6.17) 

(6.18) 


Lef Ta = Ti + R. Also define to be the first time after ri such that there are at least 3KR many left 
instructions among the random walk paths C(x,ti+i) • • • C{x,t 4 )- Clearly r4 > ri + 3KR and hence T 4 > T 3 . 
Also notice that using Lemma 23, on the event that M* (x) exists, 

M*(x) G [r3,r4]. 
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We now bound the maximum number of sleepy particles with first label x ever to be on the interval 
((x — l)K, (x + between times [rs, r 4 ], i.e. 

max Sx,x+i(t) + Sx,x-i(t). 

Condition on %{ti) and consider Clearly |r?f‘^(Ti)| < 2K, and all these particles have parts of 

random walk paths attached to it, which have not been revealed yet, and has the laws of independent lazy 
random walks (started at their respective current locations), and killed on hitting {(x — l)K, (x + 1 )K]. 
Choose A small enough such that with probability at least 1 — e all these remaining random walk paths 
corresponding to the particles in r/f'^(Ti) will not take any lazy steps at all (i.e., before hitting (x — 1)K or 
(x + 1)K). Further, choose A sufficiently small so that with probability 1 — e, Cx,j is either a left instruction 
or a right instruction for all j G {ri + 1, ri + 2,..., T 4 }. Again this is possible since for A sufficiently small, 
the chance of observing either a left instruction or a right instruction can be made arbitrarily close to one. 

By choice of R with probability at least 1 — 2e there has been at least a left instruction and a right 
instruction in the interval (ti, T 3 ]. Let 

Si = min{A: > ri : (x,k is a left instruction}, S 2 = minj/c > ri : Cx,k is aright instruction} 

and set s = si V S 2 . Now consider By definition all the intermediate sleepy particles from 

have been woken up. Now since with probability 1 — e none of their remaining paths contain any sleep 
instruction. Moreover with failure probability at most e no random walk (xj for y G [si, T 4 ] will have no lazy 
steps and hence no new particles fall asleep. Hence with probability at least 1 — 4e, we have Sx,x-i{t) = 0 
for all t G [s, r 4 ]. Taking care of Sx,x+i{t) in a similar manner it follows that with conditional probability 
at least 1 — 8 e 

max Sx,x-i{t) + Sx,x+i{t) = 0 . (6.19) 

T3<t<T4 

With this preparation, let us now return to the proof of (6.15). 


Proof of (6.15). Fix j > 0. We want to prove 


E 


AS: 


j>3 + ^~^Sj^j — l 


nPnAl] \^j-i 




Recall from previous discussion and using (6.9) ij is measurable. Recall also the filtration ^j{t) from 
(6.17) and the stopping time ri = Ti(j) from (6.18). Define ’■= A^j-i U (^j(ri). We compufe 


E 


ASj,j+i+Sj,j-i] 1 j) n D n A] 




Lef Bi denofe fhe evenf fhaf M^{j) exisfs and Sjj-i{M^{j)) + Sjj+i(M*(j)) > 0. If follows from (6.19) 
fhaf for 82 := R-iv, j) CiV Ci ACi Bi we have 

F[B2 I < 8e. (6.20) 

Now sterling from we fry fo bound fhe probabilily fhaf R{y,j) holds, i.e., in particular, 

ej{M,{j)+yj)=ij{M4j)). 

Clearly for fhis fo happen, cannof be a lefl insfrucfion for any k G [M*(j) + 1, M^{j) + yf - Lef B^ 
denofe fhis evenf. Clearly B^, is independenf of Bi and and for A sufficienlly small, 

¥[ 83 ] < {0.5l)yp (6.21) 


Now consider a consfanf C fo be specified laler (whose purpose will be clear soon). Lef B 4 denofe fhe 
even! fhaf none of fhe random walks C(j,M*(j)+i)) • • • > C{j,M*{j)+c) takes a lazy step (see (4.1)). By faking 
A = A(C') small enough 


¥[84 I — ^ £• 
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( 6 . 22 ) 






We then conclude, 


E 


eKj+i+Sjj-ili j)nvr\A] I 


( e2'^(0.51)% ifyj>C 

\ 9ee^^ + (1 - 8e) (0.51)2^^ if yj < C 


(6.23) 


The first case in the above computation is straightforward: we use the fact + Sjj-i is at most 2K 

and from (6.21), F[TZ{y;j) n 2? n ^ || < ( 0 . 51 )^t. For the second case, notice for yj < C, 


E 


e[4j+i+s,.,-i]i j)c\Vr\A] I < lP(i32 U Bi \ + F{B^ n 




The second conclusion in (6.23) now follows from (6.20), (6.21) and (6.22). Now e in (6.23) can be chosen 
to be small enough (depending on C) such that, averaging over i^j (Ti) we get that 


E 


e['^^■’^■+l+^^’^-l]l[7^(y; j) nPn^] \ ^j-i] < (0.52)^T 


This completes the proof of (6.15) and (6.16) can be established similarly. 
Now we complete the proofs of Lemma 26 and Lemma 22. 


□ 


Proof of Lemma 26. Using (6.15) and (6.16) and the naive bound of 2K on 5o,i + 5*0,-i as discussed right 
after (6.16), we have 

Thus using (6.14) and Markov’s inequality we get, 

r\V r\ Ar) < 

completing the proof of Lemma 26. □ 


Proof of Lemma 22. Using Lemma 26 and summing over all possible sequences y (see (6.7)) and using 
Observation 25 we get 

,^,,^(P(ao,/o+,/o-) n A) < ^ ^^.24) 

for some c > 0 choosing K = iT(/u) to be large enough, where the last inequality holds for large enough r. 
This establishes (6.4) and completes the proof. □ 

Next we analyse the remaining case i.e. when M(0) is large. 

6.2. M(0) is large. Recall the renormalized odometer M(-) from (5.1). To complete the proof of Lemma 
18, we still need to deal with the case when M (0) > r®. We have the following lemma. 

Lemma 27. In the set up of Lemma 18, we have 

^*,^(M(0) > r®, A) < (6.25) 

for some constant c > 0. 


Proof. For z G Z and T > 1, let L'-{T) (resp. R'fT)) denote the number of random walks among £ 

{!,... ,r}} that reach (i — 1)K (resp. (z + l)iT). Clearly, for each i and T, L[{T) ~ Bin(r, ^). By 
Azuma-Hoeffding inequality and a union bound, there exists a constant c > 0 such that with probability at 
least 1 — e~^^, the following event occurs. 


n := 



—2r 2r'l 

IT’’"’Aj 


vr > r^|L'(r) 


R!,{T)\ < 



(6.26) 
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We show below that onT-L r\ A, M (0) < r®, clearly this implies the statement of the lemma. 
Observe that on n {M(0) > r®}, we have 


f ?(0 ^ 1 ) > 


M(0) 

2 


M(0)2/3 

2 


- 2K. 


Also recall that, by the assumption on i), < 2^r, \F^ \ < 2^r for each i. Using (6.2) this implies that 

for any i > 0, 


L{i i + 1) > R{i 

We can thus conclude that for each i > 0, 


i + 1) — 4/ur. 


> m(z)/2 - M{ifl^l2 -2K- Afir. 


The above simplifies to, 


M{i + 1) > M{i) - max(M(f), -2K - 4//r. 


Let z* G {1, 2,..., ^} be the index at which M(-) is maximized. It follows that M( I) > M{u) - 
M(i^Yl^.r — {2K + 4/i)r > 0 for r sufficiently large, since by hypothesis M(L ) > M(0) > r®. Hence 
T-Ln An {M{0) > r®} = 0, and we are done. □ 


Finally we put everything together to establish Lemma 18. 

6.3. Proof of Lemma 18. By the assumption on f), notice that \Fq\, |Fg“| < 2//r. Hence taking a union 
bound over all triples (oq, fo,fo) satisfying the conditions in the statement of Lemma 22, it follows from 
Lemma 22 that 

P(A,M(0) < r®) < 

Lemma 18 now follows from the above bound and Lemma 27. □ 


7. Concluding Remarks and Open Questions 

In this paper, we have shown that ARW on Z remains in active phase even starting with arbitrarily low 
density of particles provided the sleep rate is sufficiently small. In particular this implies that the critical 
density /rc(A) < 1 for small enough sleep rate A. However our understanding of the process is still far from 
complete. Investigating the following seem the natural next step. 

(1) What happens when A is large? It is believed that for any A finite, /Uc(A) < 1. However the only 
known results in this direction is that /Zc(A) < 1 for all A ( [1,26,30]) and = 1 for A = oo [6]. 
Observe that, our arguments, though wasteful at many places, uses crucially that A can be taken to 
be arbitrarily small, it is not clear whether it is possible to improve upon this. 

(2) ARW on Z'^: Does Theorem 1 hold for ARW on Z'^ for d > 21 It is believed that the answer to this 
question is affirmative. 

(3) Critical density for SSM: For the Stochastic Sandpile Model, is the critical density Hc strictly less 
than one? As mentioned before, numerical evidence suggests an affirmative answer to this question, 
while the best known rigorous bound in [26] gives |Uc < 1- A strict inequality here would be a 
substantial progress because of the same reasons as in Remark 3. 
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